Symmetry Constraints and Variational Principles in Diffusion Quantum Monte Carlo 

Calculations of Excited-State Energies 



W. M. C. Foulkes 

Blackett Laboratory, Imperial College of Science, Technology and Medicine, Prince Consort Road, London SW7 2BZ, England 

Randolph Q. HoodO and R. J. Needs 
Cavendish Laboratory, Madingley Road, Cambridge CB3 OHE, England 
(February 1, 2008) 

Fixed-node diffusion Monte Carfo (DMC) is a stochastic algorithm for finding the lowest energy many-fermion 
wave function with the same nodal surface as a chosen trial function. It has proved itself among the most accurate 
methods available for calculating many-electron ground states, and is one of the few approaches that can be applied 
to systems large enough to act as realistic models of solids. In attempts to use fixed-node DMC for excited-state 
calculations, it has often been assumed that the DMC energy must be greater than or equal to the energy of 
the lowest exact eigenfunction with the same symmetry as the trial function. We show that this assumption is 
not justified unless the trial function transforms according to a one-dimensional irreducible representation of the 
symmetry group of the Hamiltonian. If the trial function transforms according to a multi-dimensional irreducible 
representation, corresponding to a degenerate energy level, the DMC energy may lie below the energy of the 
lowest eigenstate of that symmetry. Weaker variational bounds may then be obtained by choosing trial functions 
transforming according to one- dimensional irreducible representations of subgroups of the full symmetry group. 



PACS: 71.10.-w, 71.15.Th, 02.70.Lq, 31.25.-v 
I. INTRODUCTION 

Quantum Monte Carlo (QMC) methods are powerful 
and general tools for calculating the ground-state elec- 
tronic properties of atoms, molecules, and solids. Since 
the computational cost increases only as the cube of the 
number of particles, it is possible to study systems con- 
taining hundreds of electrons subject to periodic bound- 
ary conditions. This is enough to model real condensed 
matter with surprising precision, as shown by the accu- 
racy of 0.1 eV per atom or better achieved in QMC calcu- 
lations of the cohesive energies of solids. By comparison, 
the errors in local density functional calculations of co- 
hesive energies are often of the order of 1 eV per atom. 

The two most widely used QMC methods are varia- 
tional Mwte Carlo (VMC) and diffusion Monte Carlo 
(DMC).M In VMC a trial many-electron wave function is 
chosen and expectation values are evaluated using Monte 
Carlo integration, which is more efficient than grid-based 
quadrature methods for high-dimensional integrals. Most 
VMC simulations of solids use trial wave functions con- 
taining a number of adjustable parameters, the values 
of which are determined by minimizing the energy or its 
variance. 

DMC is a stochastic method for evolving a solu- 
tion of the imaginary-time Schrddinger equation. The 
imaginary-time evolution gradually enhances the ground- 
state component of the solution relative to the excited- 
state components, but the algorithm does not maintain 
the fermionic symmetry of the starting state. The solu- 
tion therefore converges towards the overall ground state, 
which is bosonic. This difficulty is known as the sign 
problem. 



Although several exact solutions to the sign problem 
have been proposed, none has the statistical efficiency re- 
quired to study the large systems of interest to condensed 
matter physicists. Most DMC simulations therefore use 
the approximate fixed-node method,B which is numeri- 
cally stable and often very accurate. The details will be 
described in Sec. |fi|, but the basic idea is quite simple. 
A real trial many-electron wave function is chosen and 
used to define a trial nodal surface, which is the sur- 
face on which the trial function is zero and across which 
it changes sign. In a three-dimensional system contain- 
ing TV electrons, the trial wave function is a function of 
3N variables, and the trial nodal surface is 3N — 1 di- 
mensional in general. The fixed-node DMC algorithm 
maintains the nodal surface of the trial wave function, so 
enforcing the fermionic symmetry and producing the low- 
est energy many-electron wave function consistent with 
that nodal surface. 

Although VMC and DMC are principally ground-state 
methods, they can also provide some information about 
excited states. In particular, they can be used to study 
the lowest energy state of each distinct symmetry. In 
VMC this is done by choosing a trial wave function which 
possesses the required symmetry for all values of the vari- 
ational parameters. The energy obtained after optimiz- 
ing the trial function is therefore greater than or equal 
to the eigenvalue of the lowest energy eigenstate of that 
symmetry. A similar technique is also used in DMC, al- 
though this is much harder to justify. The problem is 
that the DMC trial function is only used to define the 
trial nodal surface, which may not be sufficient to fix 
the symmetry of the state produced by the stochastic 
DMC algorithm. In any case, practical tests have shown 
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that this approach often gives excellent results. Exam- 
ples are the stuctaof excitations of the hydrogen molecule 
by Grimes eJ-al.p and calculations of excitation energies 
in diamondM and silicon.Q 

If the trial function used in an excited-state DMC sim- 
ulation has no definite symmetry, the only certainty is 
that the DMC energy must be greater than or equal to 
the many-electron ground-state energy.EI In cases when 
the DMC trial function does have a definite symmetry, 
however, it is normally assumed that the fixed-node DMC 
solution has the same symmetry as the trial function, 
and hence that the DMC energy is greater than or equal 
to the eigenvalue of the lowest energy eigenstate of that 
symmetry. This symmetry-constrained variational prin- 
ciple is widely accepted, but we show by constructing a 
specific example that it is not always correct: the fixed- 
node DMC solution need not have the same symmetry as 
the trial function; and the fixed-node DMC energy may 
be lower than the energy of the lowest exact eigenstate 
of that symmetry. 

The symmetry-constrained DMC variational principle 
is guaranteed to hold only when the trial function trans- 
forms according to a one-dimensional irreducible repre- 
sentation of the symmetry group of the Hamiltonian. The 
corresponding eigenstate is then non-degenerate, or has 
only accidental degeneracies. If the trial function trans- 
forms according to a multi-dimensional irreducible rep- 
resentation, corresponding to a degenerate energy level, 
the DMC energy may lie below the energy of the low- 
est eigenstate of that symmetry. In such cases a weaker 
variational principle may be obtained by choosing a trial 
function that transforms according to a one-dimensional 
irreducible representation of a subgroup of the full sym- 
metry group. The DMC energy is then greater than or 
equal to the eigenvalue of the lowest exact eigenstate with 
that subgroup symmetry. This provides a strict vari- 
ational lower bound for the DMC energy, but one that 
usually lies below the energy of the degenerate eigenstate 
of interest. 

As an example, consider the case of a crystalline solid. 
Any trial function with a definite crystal momentum k 
satisfies the many-electron version of Bloch's theorem 
and so transforms according to a one-dimensional irre- 
ducible representation of the translation group, which is 
a subgroup of the full symmetry group. The weaker vari- 
ational principle therefore guarantees that the DMC en- 
ergy must be greater than or equal to the energy of the 
lowest exact eigenstate with crystal momentum k. Un- 
fortunately, most Bloch states are complex and so cannot 
be used as fixed-node DMC trial functions. Real linear 
combinations of Bloch functions and their complex conju- 
gates can be used instead, but in most cases these do not 
transform according to one-dimensional irreducible rep- 
resentations and do not lead to useful variational princi- 
ples. This is illustrated in Sec. VIII, where we show that 



The weaker variational principle is useful, but relies on 
a very careful choice of trial functions and cannot explain 
all the past successes of the fixed-node DMC method for 
excited states. The real explanation of these successes, 
we believe, is that although the DMC algorithm does not 
always preserve the symmetry of the trial function, the 
imposed nodal surface acts as such a strong restriction 
that the DMC solution cannot stray "too far" from that 
symmetry. The calculated energy is therefore close to 
the variational value that would have been obtained if 
the symmetry had been preserved. In cases when the 
excited state of interest satisfies the strong variational 
principle, the errors in the ground and excited state en- 
ergies are guaranteed to have the same sign and tend to 
cancel, so improving the accuracy of the calculated en- 
ergy difference. 

The rest of this paper is organized as follows: Sec. || 
contains brief explanations of the DMC method and 
the fixed-node approximation for ground states; Sec. |l| 
shows that no general symmetry-constrained variational 
principle exists; Sec. IV shows that a variational theorem 
holds for the lowest energy state of each symmetry pro- 
vided that the trial function transforms according to a 
one-dimensional irreducible representation of the group 
of spatial transformations of the Hamiltonian; Sec. [y| in- 
troduces a weaker variational principle which may give 
energy bounds even when the trial function transforms 
according to an irreducible representation of dimension 
greater than one; Sec. VI shows that a generalization 



of the tiling theoremEl holds in every case when we can 
prove tha t th e DM C energy obeys a variational bound; 
Sees. VII and VIII give examples which demonstrate the 
absence of a general symmetry-constrained variational 
principle and illustrate the application of the weaker vari- 
ational principle; and Sec. IX summarizes and concludes. 



II. FIXED-NODE DIFFUSION MONTE CARLO 
FOR GROUND STATES 

In this section we summarize the known results con- 
cerning the application of the fixed-node DMC method 
to ground states Etu The aim is to evaluate expectation 
values with an antisymmetric wave function &(X), where 
X = (x\, #2) • • • ,xn) lists the coordinates of all N elec- 
trons, and Xi — (ri,s Z i) specifies the position and spin 
projection of electron i. We choose wave functions with 
a fixed total S z — Yli=i s zi- The expectation value of a 
spin-independent symmetric operator A(R) is given by 



(A) = 



Y^sf $*(X)i(R)$(X)dR 
Y /S J^(X)<i>(X)dR 



(1) 



the DMC energy obtained using such a trial function may 
lie below the energy of the lowest eigenstate with crystal 
momentum k. 



where R = (ri, r2, . . . , tn). For each spin configu- 
ration, S, the antisymmetric wave function $(X) = 
&{xi, X2, ■ ■ ■ , %n) may be replaced by a version with per- 
muted arguments, $(xi x , Xi 2 , . . . , Xi N ), where the permu- 
tation is chosen such that the first arguments are 
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spin-up and the last N± = N — N-\ are spin-down. Since 
R is a dummy variable, we can relabel (r^ , r i2 , . . . , r iN ) 
as (ri,r2, . . . , rjy), after which the sums over spin con- 
figurations can be removed. The expectation value may 
then be written as 

= f$«(R)i(R)$(R)dR 
J$*(R)$(R)dR ' 

where 

$(ri,r 2 ,...,rjv) = 
$(ri, t; r 2 ,T; • • ■ ; rjv, , T; rjv T +i, 4; • ■ ■ ; rjv, I) • (3) 

The function $(ri, . . . , rjv) is antisymmetric under inter- 
change of any two of the up-spin arguments ri , . . . , r jv t or 
any two of the down-spin arguments + i, . . . , rjv T +jv.) 
but has no definite symmetry under interchange of up- 
and down-spin arguments. It still obeys the Schrodinger 
equation, but the up- and down-spin electrons are now 
treated as distinguishable, allowing us to avoid explicit 
reference to the spin variables. For simplicity, all wave 
functions in the following discussion will be chosen to be 
of this type. 

DMC is a method for solving the imaginary-time 
Schrodinger equation, 

where is shorthand for ^i=i ^1- The potential V^(R) 
includes the electron-electron interactions as well as the 
external potential terms, and is assumed to be a local 
function of R. It has Coulomb singularities whenever 
two charged particles approach each other, but is finite 
everywhere else. The constant energy shift Es has been 
introduced to set a convenient zero of energy as explained 
below. The variable r (which is real) is usually called the 
imaginary time, but we will often abbreviate "imaginary 
time" to "time" in what follows. 

If the starting state ^(R, r = 0) is written as a linear 
combination of energy eigenfunctions, 

*(R,r = 0) = ^c^(R) , (5) 

i 

the large r limit of the solution of Eq. (^) takes the form, 

*(R,r->oo) =c^(R)e~ (Bi ^ s)r , (6) 

where ^(R) is the lowest energy eigenfunction appear- 
ing in Eq. (|B|) and Ei is the corresponding eigenvalue. 
In applications to many-electron systems, ^(R, r = 0) 
is antisymmetric and ^(R) is usually the many-electron 
ground state ^o(R)- We will assume that the Hamil- 
tonian H is real (possesses time- reversal symmetry), in 
which case ^(R) may also be chosen real and may be ob- 
tained by following the imaginary-time evolution of a real 



function ^(R, r). If the Hamiltonian is complex (does 
not possess time- reversal symmetry), as when there is an 
applied magnetic field, the fixed-node DMC method dis- 
cussed in this paper does not apply and it 4s necessary to 
use the fixed-phase method of Ortiz et alB 

DMC solves Eq. (||) using a stochastic algorithm, the 
efficiency of which is much, improved by an importance 
sampling transformation.EHI A real trial wave function 
$t(R) is chosen and Eq. (Q) is recast in terms of the 
product 

/(R,r) = *(R,r)$ T (R) (7) 

which is also real. After some straightforward algebra, 
/(R, t) is found to satisfy the equation, 

- |l/(R,r) = - iv£/(R,r) + V R • [F(R)/(R,t)] 
+ (£ L (R)-£ s )/(R,r), (8) 

where -El(R) = < I , x 1 (^')^ < ^t(R) is known as the local 
energy and F(R) = <I > 1 ; 1 (R)Vr$t(R) as the quantum 
force. If the chosen trial wave function is close to the ex- 
act ground-state wave function, the local energy is nearly 
constant and the statistical efficiency of the algorithm is 
optimized. 

If /(R, r) is constrained to be positive, Eq. (g) may be 
interpreted as describing the time evolution of the density 
of a population of "random walkers" multiplying or dying 
out as they diffuse and drift through a 3-ZV-dimensional 
"configuration space" . The constraint that J is positive 
is known as the fixed-node approximation^ because it 
forces the nodal surface of ^(R, t) to be the same as that 
of <3?t(R)- In practice, the distribution / is represented 
by a few hundred walkers propagating stochastically ac- 
cording to rules derived from Eq. (pj). The growth or 
decay rate of the total number of walkers in the simula- 
tion depends on the average value of -El(R) — Eg, and 
the constant energy shift E$ is chosen to ensure that the 
population remains stable on average. 

The initial walker positions are normally picked from 
the probability distribution /(R, r = 0) = (<&t(R)) 2 , 
resulting in an initial population scattered throughout 
the configuration space. The nodal surface of <&t(R) 
divides this space into different nodal pockets, among 
which the walkers are distributed. During the simulation, 
the walkers never cross the fixed nodal surface separating 
one pocket from another, and so the fixed-node DMC 
algorithm proceeds independently in each pocket. 

In the long-time limit, the probability density / of the 
walkers within nodal pocket v a becomes proportional to 
Q ,(R)<i> T (R), where the "pocket ground state" Q (R) is 
the lowest energy real normalized wave function which 
is zero outside v a and satisfies the fixed-node boundary 
conditions on the surface of v a . This function generally 
has gradient discontinuities across the surface of v a , and 
the action of the kinetic energy operator on these dis- 
continuities produces delta function terms which will be 
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denoted S a . The pocket ground state therefore satisfies 
the equations: 



H</> a (R) =e Q Q (R)+5 o 
<^ q (R)$t(R) > 

K (R) = 



when R G v a 

when R ^ w Q . (9) 



Since the walker density is positive, the sign of <j) a (R) 
within v a is the same as that of $t(R-), implying that 
4> a (R) and <I>t(R) have a non-zero overlap. The form 
of <j) a (R) does not depend on the details of <&t(R), but 
only on the shape of its nodal surface. The energy given 
by the DMC simulation in nodal pocket v a is equal to 
the pocket eigenvalue e a . This need not equal the exact 
ground-state eigenvalue Eq unless the trial nodal surface 
is the same as that of the exact ground state 'I'o(R-), i n 
which case 4> a (R) is proportional to 'I'o(R) within. u a . 

These results r-were used by Reynolds et alS3 and 
Moskowitz et alS3 to prove that the fixed-node DMC 
energy is greater than or equal to the exact ground-state 
energy. Reynolds' proof starts from the solution in a sin- 
gle nodal pocket, and uses the permutations P that do 
not flip spins to construct a real wave function antisym- 
metric with respect to interchanges of electron coordi- 
nates of the same spin, 

*a(R) = t^ttt ]T(-1) P Q (PR) = MctQl) • (10) 
T ■ J- ■ p 

The function ^(R) cannot equal zero everywhere since, 
at any point R, all terms contributing to the sum in 
Eq. ( |l0| ) have the same sign. This follows from the anti- 
symmetry of the trial function $t (R-) , which guarantees 
that all permutations P for which the point PR lies in 
nodal pocket v a have the same parity (positive if the 
sign of the trial function at R is the same as the sign of 
the trial function in v a , negative otherwise). The func- 
tion 4' ct (R) is however equal to zero everywhere on the 
trial nodal surface. This can be deduced from the knowl- 
edge that 4> a (R) is equal to zero there (as well as in 
many other places), and that the nodal surfaces of an- 
tisymmetric functions such as $t(R) are not altered by 
permutations. 

The real antisymmetric function ^ Q (R) is now substi- 
tuted into the standard quantum mechanical variational 
principle to give, 



Eq < 



J y a Hy a dR _ J m a HA(t) a dR 



J * Q * a dR 

J ^ a H<j) a dR 
J V a (j> a dR 



J ^ a A(f) a dR 



(11) 



where we have used the fact that A commutes with if, 
that it is self-adjoint, and that it is idempotent (so that 
A^ a (R) = ^q(R))- The delta function terms appear- 
ing in H <j) a do not contribute to the energy expectation 
value because they occur on the fixed nodal surface where 



^(R) = 0. If the nodal surface of $t(R) is the same as 
the nodal surface of the exact ground state, the equality 
holds and the pocket eigenvalue e a is equal to E ; but 
if <&t(R) does not have the correct nodal surface then 
e a > Eq. The energy e a produced by the DMC simula- 
tion in nodal pocket v a is therefore minimized and equal 
to the exact ground-state energy when the trial nodal 
surface is exact. Since e a is also expected to depend 
smoothly on the shape of the nodal surface, it follows 
that the error in e a is in general second order in the er- 
ror in the nodal surface. 

For systems with reasonably well behaved local poten- 
tials and real ground-state wave functions, the tiling the- 
oremB states that all the ground-state nodal pockets are 
related by permutation symmetry. A derivation of this 
theorem appears in Sec. Wl. The tiling theorem holds 
even when the ground state is degenerate, in which case 
every possible real linear combination of the degenerate 
ground states possesses the tiling property. Many DMC 
simulations use trial wave functions with the same nodal 
surface as the density-functional ground state. Since 
the density-functional Hamiltonian has a local potential, 
such trial states always satisfy the tiling theorem. This 
guarantees that the value of e a is the same in every nodal 
pocket, and hence that the energy obtained in the fixed- 
node DMC simulation cannot depend on how the walkers 
arc distributed among the pockets. 

In cases when the trial function does not possess the 
tiling property, the walker population grows most rapidly 
in nodal pockets with low values of e a (or, equivalently, 
low average values of £l(R) — -Es)- Since E$ is cho- 
sen such that the total walker population stays roughly 
constant, and since fixed- node DMC walkers never cross 
from one nodal pocket into another, the walkers in less 
favorable pockets are gradually annihilated. Although 
many different nodal pockets may have contained walkers 
initially, the population becomes more and more concen- 
trated in the pocket or pockets with the lowest value of 
e Q , and it becomes more and more likely that only these 
pockets will be occupied. The DMC energy is therefore 
almost certain to converge to the lowest of the pocket 
eigenvalues of the nodal volumes initially occupied with 
walkers, £"dmc — rnin a (e Q ). This means that the energy 
obtained in a DMC simulation which does not possess the 
tiling property may depend on the initial distribution of 
walkers. 



III. ABSENCE OF A GENERAL VARIATIONAL 
PRINCIPLE FOR FIXED-NODE DMC 
CALCULATIONS OF EXCITED STATES 

In a few cases, the exact nodal surface of an excited- 
state wave function can be determined using symmetry 
arguments alone. An example is the first excited state of 
a particle in a one-dimensional square well, which has a 
single nodal point located at the well center. The trial 
wave function can then be chosen to have exactly the 
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same nodal surface as the excited state, and the fixed- 
node DMC simulation gives the exact excited-state eigen- 
value. This result holds whether or not the trial wave 
function is orthogonal to all the lower energy eigenstates. 

In practice, however, the exact nodal surface is rarely 
determined by symmetry alone, and it is rarely possible 
to choose a trial function with exactly the same nodal 
surface as the excited state of interest. Furthermore, the 
wave function of an arbitrary excited state need not pos- 
sess the tiling property, and so the DMC energy may 
depend on the initial distribution of the walkers. 

As an example, consider a hydrogen atom in its 2s 
state, ^2s(r), with eigenvalue i?2s- The exact nodal sur- 
face is a sphere of radius ro, the value of which cannot be 
determined using symmetry arguments. If the imposed 
nodal surface is exact, the pocket eigenvalues in the in- 
ner and outer pockets will both be exactly equal to E2 S ', 
but if the nodal surface is a sphere of radius a ^ ro, 
the pocket eigenvalue e> of the wave function , 5 > in the 
outer pocket will not equal the pocket eigenvalue e< of 
the wave function \& < in the inner pocket. 

Consider the case when the fixed node is too close to 
the nucleus, a < r$. According to the variational prin- 
ciple applied to the outer pocket, the value of e> must 
be bounded above by the energy expectation value of the 
trial function, 



*>(r) 



a < r < ?'o 

*2 S (0 r > r . 



(12) 



This expectation value is the exact 2s eigenvalue -E^s- 
Clearly, we could construct a lower energy trial function 
by removing the kink in ^^(r) at r — ro, and hence 
e> < i?2s- Similarly, we can take the pocket ground state 
from the inner pocket, and use it as a variational 
trial function within the exact 2s nodal surface r = ro, 



$<(r) 



*<(r) r < a 

a < r < r . 



(13) 



In this case the energy expectation value of the trial func- 
tion is e< , while the minimum possible expectation value 
for a wave function with a node at ro is E2 S - Again, we 
could construct a lower energy trial function by removing 
the kink at r = a, and hence e< > E% s . Note that varia- 
tional arguments like these can be used whenever the ex- 
act nodal pocket completely encloses the trial pocket or 
vice-versa, irrespective of dim ension or symmetry. This 
will prove useful in Sec. VIII . 

The last paragraph showed that if a < ro then e> < 
E2 S < £<; if a > fQi & similar derivation gives e< < 
E2s < £>• The two pockets have different energies un- 
less a = ro. How does this affect the fixed- node DMC 
algorithm? As long as the lower energy pocket contains 
plenty of walkers initially, the walker population in the 
higher energy pocket will almost certainly die out. The 
fixed-node DMC energy will then tend to the pocket 
eigenvalue of the lower energy pocket, which is always 
less than or equal to E2 S - The DMC energy is therefore 



maximized when a — ro. If a is increased through ro, 
the r — ► oo DMC walker population switches from the 
outer nodal pocket to the inner one, and the slope of the 
graph of the DMC energy versus a changes discontinu- 
ously. It follows that the error in the fixed-node DMC 
energy is first order in the error in the nodal surface, not 
second order as it is for the ground state. This example 
shows that there is no variational principle when fixed- 
node DMC is used to study general excited states: the 
error in the DMC energy of the excited state may in- 
crease linearly with the error in the nodal surface; and 
the DMC energy need not be minimized when the nodal 
surface is exact. 



IV. THE LOWEST ENERGY EIGENSTATE OF 
EACH SYMMETRY 

We now address the question of whether there is a vari- 
ational principle for fixed-node DMC simulations of the 
lowest energy state of each symmetry. We denote by Q 
the group of spatial transformations T (combinations of 
rotations, reflections, translations, and inversions of all 
electrons simultaneously) which leave the many-electron 
Hamiltonian invariant. For simplicity, we consider this 
group to be finitaia and of order g. The full symme- 
try group of the Hamiltonian is the direct product of 
Q with the permutation and time- reversal groups. All 
the symmetry arguments in this paper can easily be re- 
cast in terms of the full symmetry group instead of Q 1 
but real many-electron wave functions have such simple 
time-reversal and permutation symmetries (they trans- 
form according to the one-dimensional identity represen- 
tation of the time-reversal group and the one-dimensional 
antisymmetric representation of the permutation group) 
that this is unnecessary. Since the arguments based on 
the spatial group Q are somewhat easier to grasp, we 
choose to work with this group in what follows. 

We begin by showing that a symmetry-constrained 
variational principle holds whenever the chosen trial 
function is real and has an invariant nodal surface. This 
theorem was stated without proof in a paper by Caffarel 
and Claverie.EJ Here we sketch a proof. 

A real trial function $t (R) is said to have an invariant 
nodal surface if the transformed function, 



Q(T)$ T (R) = $ T {T- l R) , 



(14) 



has the same nodal surface as $t(R) for all coordinate 
transformations T 6 Q. Note that the nodal surface is 
defined as the surface on which $t (R) is zero and across 
which it changes sign; <!>t(R) may also have additional 
zeros where the sign does not change, but these are not 
on the nodal surface and need not be invariant. 

The proof relies on the properties of the function x(R) 
defined by, 
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on the trial nodal surface 
x(R) = < +1 in nodal pockets where <&t > 
— 1 in nodal pockets where <&t < 



(15) 



Given any point R not on the nodal surface, and any 
spatial transformation T S G, it is clear that x(R) = 
77X(T _1 R), where r\ — ±1. Furthermore, since the func- 
tions x(R) an d x(^" _1 R) have the same nodal surface 
and so change sign together as R changes, the sign of rj 
must be independent of R. This shows that all symme- 
tries T £ Q either leave x(R) unchanged or multiply it 
by —1, from which it follows that x(R) transforms ac- 
cording to a one-dimensional irreducible representation 

t t q of g. 

We can now adapt the proof of the ground-state varia- 
tional principle given in Sec. [n| Take the pocket ground 
state 4> a and antisymmetrize it to obtain a function $ a as 
in Eq. pQl . Now apply the group theoretical projection 
operatorT-3 



(16) 



Teg 



where Tg(T) is the one-by-one matrix representing T. 
The application of V r produces a new antisymmet- 
ric state ^! T a which transforms according to the one- 
dimensional irreducible representation Tg. Since the 
original pocket ground state <p a (R) had a non-zero over- 
lap with x(R)) which is itself an antisymmetric function 
of symmetry Tg, the antisymmetric Tg component pro- 
jected out by applying A and then V r to <j) a cannot be 
zero. 

The energy expectation value of ^> r a is greater than or 
equal to the energy Eq of the lowest antisymmetric state 
of symmetry T r g, 



El < 



J V r a *M r a dR _ J W*HV r A$ a dR 



J dR 

J V r a *H0 a dR 
J ^ r a *(t) a dR 



J $> r a *V r A<t> a dR 



(17) 



In analogy with the ground-state proof of Sec. 0, these 
manipulations rely on the self-adjointness and idempo- 
tency of the operators A and V r , both of which commute 
with H. It is also important that (R) is equal to zero 
everywhere on the surface of v a , so that the delta func- 
tion terms in the expression for H cf> a (R) from Eq. (^|) do 
not contribute to the expectation value. This is guaran- 
teed because the fixed nodal surface is invariant. 

Note that the symmetry of the trial function played 
no part in this derivation; the only thing that mattered 
was the invariance of the trial nodal surface. In most 
DMC simulations, however, the trial function does have a 
definite symmetry, and so transforms according to a spe- 
cific irreducible representation of Q. Given that functions 
transforming according to different irreducible represen- 
tations are orthogonal, and that the definition of x(R) 



ensures that (xI'&t) > 0, it follows that any trial function 
with a definite symmetry and an invariant nodal surface 
must transform according to the same one-dimensional ir- 
reducible representation as the corresponding x(R)- This 
implies that a real trial function transforming according 
to an irreducible representation of dimension greater than 
one cannot have an invariant nodal surface. 

The converse is also true: any real trial function <1>t (R) 
transforming according to a one-dimensional irreducible 
representation must have an invariant nodal surface. 
This follows from the transformation law, 



Q(T)$ T (R) = r£(T)$ T (R) , 



(18) 



and the observation that the real normalized function 
$t(R) remains real and normalized under all transfor- 
mations in Q. The one- by-one matrix Tg(T) is therefore 
equal to ±1, and the nodal surface of <5(T)$t(R) is the 
same as that of <1>t(R)- 

Putting everything together, we can now conclude that 
whenever the real trial function transforms according to 
a one-dimensional irreducible representation of Q, the 
nodal surface is invariant, and the DMC energy is greater 
than or equal to the energy of the lowest exact eigen- 
function with that symmetry. This is the symmetry- 
constrained variational principle mentioned in the intro- 
duction. If the trial function transforms according to 
an irreducible representation of dimension greater than 
one, the nodal surface cannot be invariant, and the delta 
functions produced when H is applied to <j) a need not all 
occur where 4'™ = 0. The above proof of the symmetry- 
constrained variational principle therefore breaks down. 
In Sec. VIII we give an example of a system with a trial 
function transforming according to an irreducible repre- 
sentation Tg of dimension d r = 2 for which -Edmc < E 7 Q . 

The restriction of the variational principle to trial 
functions transforming according to one-dimensional ir- 
reducible representations can be understood in a very 
simple way. The problem being solved in a fixed-node 
DMC simulation is not the imaginary-time Schrodinger 
equation of Eq. (Q) , but the imaginary-time Schrodinger 
equation subject to the additional boundary conditions 
specified by the trial nodal surface. As a result, the rel- 
evant symmetry group is not the group Q of symmetries 
of H, but the subgroup C/fn that leaves both the Hamil- 
tonian and the fixed nodal surface (boundary conditions) 
invariant. The fixed-node eigenfunctions need only con- 
form to the symmetries in C/fn, and their transforma- 
tion properties should be analyzed in terms of the irre- 
ducible representations of Gfn, not Q. If the real trial 
function transforms according to a one-dimensional ir- 
reducible representation of £/, the fixed nodal surface is 
invariant and the groups Q and C/fn are the same. The 
symmetry-constrained variational principle discussed in 
this section does not apply unless this is the case. 
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V. IRREDUCIBLE REPRESENTATIONS OF 
DIMENSION GREATER THAN ONE 

Can anything be said when the trial function trans- 
forms according to an irreducible representation of di- 
mension greater than one? Suppose that the real trial 
function l (R) transforms as the Ith row of an irre- 
ducible representation Tg with dimension d r > 1. We 
already know that the nodal surface of <I>^ t (R) cannot 
be invariant with respect to all the operations in G, but 
it may be invariant under a subset of those operations. 
Any such subset forms a proper subgroup, Gfn, which in 
general depends on the row index I. 

As in Sec. IV, we define the function x(R) which is 
equal to +1 in all nodal pockets where $^, ; (R) > 0, 
equal to —1 in all nodal pockets where $^,j(R) < 0, 
and equal to zero everywhere on the nodal surface of 
$t ;(R). By construction, x(R) transforms according to 
a one-dimensional irreducible representation of (Jfn (al- 
though not of Q). Since all the pocket ground states <p a 
have non-zero overlaps with x(R-)i they all have non-zero 
components of the same subgroup symmetry as x(R-)- 
We can therefore re-use the symmetry-projection argu- 
ment leading to Eq. ([l7]) to show that the fixed-node 
DMC energy is greater than or equal to the eigenvalue 
of the lowest energy exact eigenfunction with the same 
subgroup symmetry as x(R). This provides a rigorous 
variational principle, but the subgroup Gfn that leaves 
the trial nodal surface invariant is usually small and its ir- 
reducible representations provide a correspondingly lim- 
ited set of symmetry labels. The bounds obtained are 
therefore weak and this variational principle is of little 
use unless optimized as explained below. 

If the nodes of <&t { (R) were exact one could form a 
trial function from any real linear combination of differ- 
ent rows and always obtain the same DMC energy. When 
the nodes are not exact, however, the DMC method 
breaks this <i r -fold degeneracy, and the variational lower 
bound on the DMC energy depends on the precise linear 
combination of rows chosen in constructing the trial wave 
function. This freedom can be exploited to improve the 
weak variational principle. 

The d r functions <j>ip ; (R) (I = 1,2,..., d r ) are a ba- 
sis for the irreducible representation Tg of Q, and hence 
transform into linear combinations of each other under 
all coordinate transformations T E Q. This implies that 
they also transform into linear combinations of each other 
under all coordinate transformations T € Q s , where Q s 
is any proper subgroup of Q . The subset of the matrices 
T r g{T) for which T G Gs is therefore a representation of 
Q s . This representation is not in general irreducible. 

To enumerate the different possibilities, begin by de- 
termining all the proper subgroups Gs oiQ. For each such 
subgroup, consider the representation of Gs consisting of 
the matrices Tg(T) with T € Q s . This representation 
may be decomposed into its irreducible components with 
respect to Gs, 
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' ' ffi n q m Tg™ 



(19) 



where the positive integer n qi is the number of times 
the irreducible representation Tg of Gs appears, so d r — 
J2i n qid qi - The irreducible representations of Gs appear- 
ing in Eq. ( |l9| ) are said to be compatiblai-3 with Tg . Note 
that since the group containing only the identity element 
is always a subgroup of G, it is always possible to find at 
least one subgroup for which the reduction of Tg contains 
a one-dimensional irreducible representation Tg . 

A trial function transforming as the Ith row of Tg may 
contain components along all the compatible representa- 
tions of Gs, but it is always possible to construct linear 
combinations, 



^(R)=£cf^(R), 



(20) 



i=i 



transforming according to each particular compatible 
representation T'g . Every real function <&^ l (R) corre- 
sponding to a one-dimensional representation of Gs has 
an invariant nodal surface with respect to Gs ■ A DMC en- 
ergy calculated using the nodes of such a function there- 
fore satisfies the variational principle -Edmc > Eq . 

Clearly, the strength of the variational principle ob- 
tained depends on the choices of the subgroup Gs and the 
one-dimensional representation Tg . In some cases, one 
can find a subgroup Gs with a one-dimensional represen- 
tation Tg which is compatible with Tg but incompatible 

with all those irreducible representations Tg for which 
Eq < Eq. A DMC energy calculated using the trial 
function <i>^ (R) is then guaranteed to be greater than or 
equal to Eq. In general, if one knows the ordering of the 
energy levels beforehand, one can use the compatibility 
analysis to find the one-dimensional irreducible represen- 
tation Tg which gives the most stringent energy bound 
for the eigenvalue of interest. 

As an application of this symmetry-constrained vari- 
ational principle, consider a crystal with space group 
G- We wish to establish whether it is possible to use 
fixed-node DMC to obtain a variational estimate of the 
eigenvalue E (k) of the lowest energy eigenstate with 
crystal momentum k. The relevant subgroup Gs is the 
translation group, which is Abelian and has only one- 
dimensional irreducible representations. These are la- 
beled by the crystal momentum k, and so the statement 
of the symmetry-constrained variational principle is very 
straightforward: if the DMC trial function has crystal 
momentum k, the DMC energy is greater than or equal 
to Eo(k). It is worth noting that this statement holds 
whether or not the lowest energy state of crystal mo- 
mentum k is degenerate. In the degenerate case the trial 
function transforms according to a multi-dimensional ir- 
reducible representation of G, but it still has crystal mo- 
mentum k and so still transforms according to a one- 
dimensional irreducible representation of the translation 
subgroup Gs- This is sufficient to guarantee that the 
DMC energy is greater than or equal to i?o(k). 
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The symmetry-constrained DMC variational principle 
for the lowest energy state of crystal momentum k is 
unfortunately much less useful than it appears, because 
most Bloch states are complex and cannot be used as 
fixed- node DMC trial functions. Instead, the standard 
approach is to use a real linear combination of a Bloch 
function and its complex conjugate. (This is justified by 
the assumption of time-reversal invariance, which guar- 
antees that if is an eigenfunction with eigenvalue 
i?o(k) then so is ^J.) Since such trial functions contain 
components with two different wave vectors, k and — k, 
they do not transform according to a single irreducible 
representation of the translation group. The symmetry- 
constrained variational theorem is therefore inapplicable, 
and it is possi ble th at the DMC energy may lie below 
Eo(k). In Sec. VIII we show by means of a specific ex- 



ample that such calculations may indeed produce a DMC 
energy which is lower than -Eo(k). 

An important exception arises when the wave vector k 
is equal to half a reciprocal lattice vector, in which case 
k and — k are alternative labels for the same irreducible 
representation of the translation group. The linear com- 
bination of \l/k and is then a pure Bloch function, and 
the normal proof of the symmetry-constrained variational 
principle applies. As long as k equals half a reciprocal 
lattice vector, the DMC energy is greater than or equal 
to the energy of the lowest exact eigenstate with crystal 
momentum k. 



VI. GENERALIZATIONS OF THE TILING 
THEOREM 

In his paper on fcrmion nodes Jl Ceperley stated that 
the tiling theorem could be generalized to the case where 
there were other discrete symmetries present. A more 
precise statement is that, given a Hamiltonian with a 
reasonable local potential and a symmetry group Q, the 
tiling theorem applies to any real state which is the lowest 
energy eigenfunction of a symmetry Tg with dimension 
d r = 1. This can be demonstrated via a simple gener- 
alization of Ceperley's proof! of the ground-state tiling 
theorem. 

Consider the real antisymmetric state ^(R) which 
is the lowest energy eigenfunction transforming accord- 
ing to the one-dimensional irreducible representation Tg. 
This function may have many nodal pockets, but pick 
one at random and color it blue. Now apply a symmetry 
operator from the group containing the spatial symme- 
tries (rotations, translations, reflections and inversions) 
and the permutations. Since the nodal surface is invari- 
ant, this symmetry operator maps the region of space 
within the blue nodal pocket into itself or into one of the 
other nodal pockets. If the blue pocket is mapped into 
some other nodal pocket, this pocket is equivalent to the 
blue one by symmetry and is also colored blue. Repeat 
this process for every operator in the group, until all the 



pockets equivalent to the original one have been found 
and colored blue. 

There are now two possibilities: the blue regions may 
fill the entire configuration space, in which case the nodal 
pockets are all equivalent by symmetry and the tiling 
theorem holds; or there may be other inequivalent nodal 
pockets which have not yet been found. We can rule out 
the second possibility using the following argument. 

Assume that the blue nodal pockets do not fill the con- 
figuration space, and that the local potential V(R) and 
all its derivatives are finite except at the Coulomb sin- 
gularities occurring when two charged particles approach 
each other. As long as the system is not one-dimensional 
(in which case different arguments are required), this en- 
sures that almost every point on the nodal surface lies a 
finite distance away from the nearest singularity in the 
potential. The eigenfunction V&o^) may therefore be 
expanded as a power series with a finite radius of conver- 
gence about almost any point R s on the nodal surface. 
If the gradient of ^(R) is assumed to be zero over any 
finite area of the nodal surface surrounding R s , it can 
be shown that every term in this series has to be zero, 
and hence that ^(R) = everywhere within the radius 
of convergence. Since any solution of the Schrodinger 
equation can be analytically continued around the iso- 
lated singularities in the potential, this further implies 
that \&o(^-) i s zer0 everywhere. We therefore conclude 
that the gradient of \&o(R) mus t be non-zero almost ev- 
erywhere on the nodal surface. 

Now consider the trial function \& r (R), which is de- 
fined to equal 'I'o(R-) within the blue pockets and zero 
elsewhere. This trial function is antisymmetric and trans- 
forms according to the irreducible representation Tg , but 
has gradient discontinuities on the nodal surfaces sepa- 
rating the blue pockets from the rest of configuration 
space. It satisfies the Schrodinger-like equation, 



HV r (R) = E^ r (R) + 5 r 



(21) 



where the symbol 8 r denotes the delta functions produced 
by the action of the kinetic energy operator on the gra- 
dient discontinuities. The delta functions occur where 
v& r = and so do not affect the energy expectation value, 



mm 



(\Jrr|\J/r\ 



(22) 



We know, however, that a state which has gradient dis- 
continuities almost everywhere on a finite area of the 
nodal surface cannot be an eigenfunction unless the po- 
tential is infinite almost everywhere on that area. Since 
we are assuming that this is not the case, the function 
ty r must contain excited-state components of symmetry 
Tg and cannot have the same energy expectation value 
as the lowest energy state of that symmetry. This con- 
clusion contradicts Eq. (p2]), and so the assumption that 
the blue nodal pockets do not fill the configuration space 
must have been incorrect. All the nodal pockets of 
are therefore equivalent by symmetry. 
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This proves the tiling theorem for any real state which 
is the lowest energy eigenfunction of a symmetry Tg with 
dimension d r = 1. An obvious corollary is that there is 
also a tiling theorem for the lowest energy state trans- 
forming according to any one-dimensional irreducible 
representation of any subgroup Q s of Q. This statement 
is analogous to the weaker variational principle discussed 
in Sec. |V|. In every case when we have demonstrated 
the existence of a DMC variational principle, we have 
therefore also been able to demonstrate the existence of 
a tiling theorem. Our analysis has been restricted to 
the physically interesting case of a local potential which 
is finite everywhere except at Coulomb singularities, but 
our conclusions may be somewhat more general than this 
suggests. 

The familiar many-fermion ground-state tiling theo- 
rem may be viewed as a special case of the subgroup 
tiling theorem mentioned above. The permutation group 
is always a subgroup of the full symmetry group (which 
contains both spatial and permutation symmetries), and 
the many-fermion ground state is the lowest energy state 
which transforms according to the one-dimensional an- 
tisymmetric irreducible representation of that subgroup. 
The subgroup tiling theorem therefore guarantees that 
the many-fermion ground state possesses the tiling prop- 
erty. Note that the tiling property holds with respect to 
the permutation subgroup, not the full symmetry group. 
This means that it is only the elements of the permu- 
tation subgroup that need be applied to the initial blue 
pocket to find all equivalent pockets and turn the whole 
configuration space blue. 

The above derivation of the subgroup tiling theo- 
rem only applies to states that transform according to 
one-dimensional irreducible representations of the chosen 
subgroup G s - Such states may also transform according 
to multi-dimensional irreducible representations of the 
full symmetry group G, so the tiling theorem is not re- 
stricted to non-degenerate energy levels. In systems with 
degenerate many-fermion ground states, for example, any 
real linear combination of the degenerate ground states is 
antisymmetric and so possesses the tiling property with 
respect to the permutation subgroup. This result holds 
even though the nodal surface is not invariant under all 
the operators from the full symmetry group. 

The most important consequence of the generalized 
tiling theorem is exactly as in the ground-state case. It 
is common for a fixed-node DMC trial wave function to 
have the same nodal surface as an energy eigenfunction 
calculated using an approximate method such as local- 
dcnsity-functional theory. The approximate Hamilto- 
nian is chosen to have the same symmetries as the exact 
Hamiltonian, but may also have extra symmetries which 
are not relevant to the argument and may be ignored. If 
the approximate Hamiltonian has a reasonably well be- 
haved local potential, its eigenstates have the same tiling 
properties with respect to Q as the corresponding exact 
eigenstates. The lowest energy eigenstate transforming 
according to any one-dimensional irreducible represen- 



tation of Q or any subgroup of Q therefore satisfies the 
tiling theorem. This ensures that the energy produced 
by a DMC simulation using the nodal surface of such a 
state is independent of the initial walker distribution. 



VII. TIGHT-BINDING EXAMPLE 

This section and the next describe simple examples 
that show what happens when a trial wave function 
transforming according to an irreducible representation 
with d r > 1 is used to define the trial nodal surface for 
a DMC simulation. We find that the DMC energy may 
indeed be lower than the eigenvalue of the lowest ex- 
act eigenstate with the same symmetry as the trial func- 
tion. This demonstrates that the variational principle of 
Eq. ( |l7| ) does not apply when d r >1. 

In searching for a suitable example system, we found 
it convenient to impose the following restrictions: (i) the 
group of the Hamiltonian is finite; (ii) the trial wave func- 
tion and its nodal surface can be easily visualized; and 
(iii) the exact eigenstates of the Hamiltonian and the 
pocket ground states of Eq. (||) can be calculated with- 
out numerical error. These restrictions are unnecessary, 
but make the analysis much simpler. To satisfy criterion 
(ii), we must choose a system containing at most three 
electrons in one dimension or a single electron in two or 
three dimensions. The nodal surfaces of one-dimensional 
systems are not easily altered because they are almost 
entirely determined by the antisymmetry of the wave 
function, so we decided to concentrate on one-electron 
systems in two and three spatial dimensions. 

The first system we studied was an electron confined to 
a rectangular box in two or three dimensions. Although 
this system does have multi-dimensional irreducible rep- 
resentations, we failed to find an example in which the 
DMC energy calculated using a trial state of a given 
symmetry was less than the eigenvalue of the lowest en- 
ergy eigenstate of that symmetry. The second system 
we tried was more successful. We start by discussing a 
simple tight -binding realization of this system, and then 
in Sec. VIII present an alternative realization based on a 



separable solution of the Schrodinger equation. 

Consider a molecule containing one electron moving 
in the potential of three protons fixed at the corners 
of an equilateral triangle. The symmetpyj group of the 
Hamiltonian of this system is called Z? 3 hJi-3 A convenient 
way to generate DMC trial wave functions with specific 
symmetries is to solve the Schrodinger equation within 
a tight-binding approximation using a single spherically- 
symmetric atomic-like orbital £(l r l) centered on each pro- 
ton. Once a tight-binding eigenfunction has been found 
and used to define a trial nodal surface, the fixed-node 
DMC algorithm can be used to solve the original Hamil- 
tonian exactly subject to the fixed-node constraint. 

The ground state of the molecule is a non-degenerate 
nodeless function transforming according to the identity 







representation. Some of the excited-state eigenfunctions 
must be doubly degenerate, however, since the symmetry 
group £>3h has two two-dimensional irreducible represen- 
tations, one of which is called r 3 n . 
. . . . . 3h 

The tight-binding Hamiltonian has only three eigen- 

states, the nodeless ground state and a degenerate pair 
of excited states which can be written as the Bloch func- 
tions, 



r 3 

1 D, 



^3h 



(26) 



i27rj/3 



3=0 
2 



*-(r)=££(l r - d ,l) e ~ l27rj/3 



(23) 



3=0 



where do, di, and d2 are the position vectors of the three 
protons. These two Bloch functions transform into lin- 
ear combinations of each other under the operations of 
the point group (there is no need to include time-reversal 
symmetry) and form a basis for . Taking linear com- 
binations of V E'+ and one can form the real functions, 

* 1 (r)=C(|r-d 1 |)-€(|r-d 2 |), 

* a (r) = 2£(|r - d |) - £(|r - d x |) - £(|r - d 2 |) , (24) 

which form an alternative basis for the same irreducible 
representation. 

Fig. (Q) shows the nodal surfaces of \&i and '5 2 along 
with the expansion coefficients from Eq. (p4j). The 
nodal surface of *Sfi does not depend on the particular 
spherically-symmetric tight-binding basis functions cho- 
sen and turns out to be an exact excited-state nodal sur- 
face. The nodal surface of $2 is not exact, however, and 
its precise shape depends on the details of the atomic- 
like orbitals used in the tight-binding model. Different 
choices of £(|r|) give trial functions with the same sym- 
metries but different nodal surfaces. By changing £(|r|), 
it is possible to change the relative sizes of the two nodal 
pockets of ^2, causing one to grow at the expense of 
the other. As in the example of the 2s state of the H 
atom discussed earlier, this suggests that in many cases 
the pocket eigenvalue of one of the two nodal pockets 
will be too high while that of the other is too low. The 
DMC energy will then lie below the energy of the exact 
eigenstate of interest. In Sec. VIIl we numerically solve 



a specific example with the same symmetry properties as 
this system and observe exactly this behavior. 

Here we apply the weaker variational principle de- 
scribed in the previous section to the degenerate excited 
state with symmetry r 3 -,^ of the group Dy^. The largest 
subgroups and D3 each have six elements. If the sub- 
group -D3 is used, the compatibility relation of Eq. jTg| ) 
becomes, 

(25) 



1 D 3h - 1 D 3 



where 



is two dimensional. This result is not use- 
ful since we are seeking a reduction that contains one- 
dimensional representations. Using the subgroup C^^ re- 
sults in two one-dimensional representations 



However, any trial wave function that transforms as ei- 



ther T 5 r 



or ] is complex and therefore unsuitable for 



fixed-node DMC. 

If we use the smaller subgroup, C?2 V j which has four 
elements, we obtain the compatibility relation 



r 3 



c 2v 



C 2 v 



(27) 



from which we can construct real trial wave functions and 
apply the weaker variational principle. In fact, the trial 
functions and ^2 already have the correct transfor- 
mation properties: transforms as the one-dimensional 
representation ^c 2v ' wnuc ^2 transforms as the one- 
dimensional representation T^ 2 . The group C'2v there- 
fore preserves the nodal surfaces of \&i and Because 
the symmetry corresponding to the representation rj> 
is compatible with the ground-state symmetry as well 
as with r|, , a DMC simulation using the trial func- 
tion ^ 2 satisfies only the weaker variational principle 
£dmc > Eq, where Eq is the overall ground-state en- 
ergy. The representation T^ 2v is not compatible with 
the ground-state symmetry, however, and so a simulation 
using the trial function gives a stronger variational 
principle. For the example studied in the next section, it 
turns out that the T 3 ^ state of interest is the lowest en- 
ergy exact eigenstate with which the representation r<t, 
is compatible. The strong variational principle therefore 
applies and Sdmc > Eq, where Eq is the exact T 3 -,^ 
eigenvalue. It is important to appreciate that this is not 
a general result; it is not always possible to pick a trial 
function that maintains the strong variational principle 
for a given symmetry. 

Note that the strong variational principle applies be- 
cause the trial wave function transforms according 
to the irreducible representation of the subgroup 

C2V If we had chosen a different pair of trial functions, 
constructed by taking linear combinations of ^1 and ^2, 
both would have contained components along r^, and 
^c 2v - The only subgroup preserving the nodes would 
then have been the group of the identity, and the only 
variational principle would have been with respect to the 
overall ground-state energy, Eumc > E . This illustrates 
the general rule that the strongest variational principles 
are obtained by choosing trial functions which transform 
according to specific one-dimensional irreducible repre- 
sentations of specific subgroups of Q. 



VIII. SEPARABLE EXAMPLE 

We now present an explicit solution of a different ex- 
ample with the same symmetry group as the trian- 
gular molecule discussed above. Consider a particle of 
unit mass moving in a triangular potential in three di- 
mensions. The wave function \&(r) obeys the Schrodinger 
equation, 
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1. 



V 2 + V(r) *(r) = £#(r) 



where 



V(r) = 



cos(36>) 



(28) 



(29) 



with p, 9 and z the usual cylindrical coordinates. The 
boundary conditions are ^(p, 9, z) = for p > 1 or |z| > 
7r/2, confining the particle within a cylinder. 

Writing W(r) = i?(p)9(6')Z(z), the Schrodinger equa- 
tion separates into: 



l rf 2 Z 
2 dz 2 



S«2 



-^ + cos(30)e = .E fl e 



1 d 



dR 



P- 



E, 



E 



2p dp \ dp J \ p 
The lowest energy eigenfunction of Eq 

1 



R = ER 



Z{z) 



cos(z) 



(30) 

(31) 
(32) 

(33) 



with eigenvalue E z = 1/2. Eq. (p^) m ay be simplified 
by the substitutions r = \/2(E — E z )p and v = \/2Eg, 
which yield Bessel's equation, 



o dR dR o 9 
r 2 — j + r— + (r 2 - v 2 )R = . 
dr z dr 



(34) 



The general solutions are the Bessel and Neumann func- 
tions, J v {r) and N v (r), but only J v {r) is well behaved at 
the origin. Hence 



(35) 



with the energy E determined by the boundary condition 
at p = 1, 



JvmiV^E - E z )) = 



(36) 



Eq. (|3l| ) can be transformed into Mathieu's equationlij by 
a simple change of variables, but here we solve it numer- 
ically by expanding the eigenfunctions in (normalized) 
sines and cosines, 



6(0) = a - 



1 



/2tt 



M 

E 



oobc^) , 6n !M^r, j (37) 



and diagonalizing the corresponding Hamiltonian ma- 
trix. The results converge rapidly with M, and choosing 
M = 50 gives very accurate eigenvalues for the lowest few 
eigenstates. 

The three lowest energy angular eigenfunctions are 
shown in Fig. (0). As stated in Sec. VU, the lowest 



all elements of the symmetry group of the Hamiltonian. 
The next two states are a degenerate pair forming a basis 
for the irreducible representation r^, . One of the two 
has nodes at 9 = and 9 = ir, and will be called 



in analogy with the corresponding state from Sec. VH 



the other has nodes at 9 = ±1.7934 radians, and will 
be called ^2- Note that \&i is the lowest energy eigen- 
function transforming according to the one-dimensional 
irreducible representation of the subgroup C2 V ; it 

therefore possesses the tiling property with respect to 
C<i y . The r^ 2 symmetry of <3/2 is shared by the overall 
ground state, however, so ^2 is not the lowest eigenstate 
of that subgroup symmetry and does not satisfy a tiling 
theorem. 

Consider how DMC might be used to find the eigen- 
value of the lowest energy doublet. We do not want 
to impose the exact nodal surface since DMC would then 
give the exact answer and we would learn little about 
variational principles, so we seek a trial function with 
a different nodal surface but the same symmetry. We 
choose to generate such a trial function by solving the 
Schrodinger equation for a different triangular potential, 



H 



1 9 
-2 V 



cos(36>) cos(6(9) 

15 r M 19 



(38) 



energy eigenstate is a nodeless function invariant under 



where p is an adjustable parameter. The first three an- 
gular eigenfunctions of this Hamiltonian when p = 5 are 
shown in Fig. (||). Since the cos(6(9) term does not change 
the symmetry group, these eigenfunctions still belong to 
the same irreducible representations. However, the nodal 
angle of ^2 has moved slightly in response to the pertur- 
bation. One of the two nodal volumes of this trial wave 
function is therefore slightly too small, while the other is 
slightly too large. 

The trial nodal pocket that is too large (small) encloses 
(is enclosed by) the corresponding exact nodal pocket. 
We can therefore apply the variational argument used 
in the discussion of the 2s state of the hydrogen atom 
to show that the pocket eigenvalue from the pocket that 
is too small must be greater than the exact eigenvalue, 
while the pocket eigenvalue from the pocket that is too 
large must be less than the exact eigenvalue. As in the 
H 2s example, therefore, the DMC energy is always < 
the exact eigenvalue. The maximum of the DMC energy, 
equal to the exact r^ 3h ground-state eigenvalue, is at- 
tained only when the nodal angle is exact; and the slope 
of the graph of DMC energy against nodal angle changes 
discontinuously at this point. The error in the DMC en- 
ergy is first order, not second order, in the error in the 
nodal angle. 

This analysis is confirmed by the results of a full cal- 
culation given in Fig. (||), which shows how the angular 
pocket eigenvalues Eg of the two pockets depend on the 
angular half width of the nodal pocket centered on 9 = 0. 
As expected, Eg is too large when the nodal pocket is too 
small and vice-versa. Fig. (|5|) shows the dependence of 
the total pocket eigenvalues E on the angular half width 
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of the nodal pocket centered on 6 = 0, confirming that 
an increase in angular half width gives rise to a decrease 
in total eigenvalue and vice-versa. The numerical results 
therefore support the conclusions of the variational argu- 
ment. 



Note, finally, that Eq. ( |31[) can be interpreted as the 
Schrodinger equation for a one-dimensional crystal with 
periodic boundary conditions, in which case and $2 
are real linear combinations of Bloch waves with equal 
and opposite crystal momenta. Seen from this viewpoint, 
the degeneracy of and 'J' 2 arises from the time-reversal 
(complex-conjugation) symmetry of the real Hamiltonian 
rather than from its spatial symmetry, but the failure 
of the symmetry-constrained variational principle is still 
apparent. This confirms the assertion made in Sec. 0: 
a DMC simulation using a real trial state constructed 
from Bloch states with equal and opposite crystal mo- 
menta may yield an energy below that of the lowest exact 
eigenstate with that crystal momentum. 



IX. CONCLUSIONS 

The main lesson to be learned from this work is that 
symmetry arguments cannot be applied to fixed-node 
DMC unless the symmetries of both the Hamiltonian and 
the nodal surface of the trial wave function are taken 
into account. The fixed-node DMC algorithm solves the 
Hamiltonian subject to the boundary conditions imposed 
by the trial nodal surface, and so the relevant symme- 
try group (Jfn contains only those symmetry operations 
that leave both the Hamiltonian and the trial nodal sur- 
face (boundary conditions) invariant. This "group of the 
fixed- node Hamiltonian" is a subgroup of the more famil- 
iar "group of the Hamiltonian", Q. 

If the trial function transforms according to a real one- 
dimensional irreducible representation r r of Q, all sym- 
metry operations in Q simply multiply the trial function 
by a real number. This does not change the nodal sur- 
face and hence Q = Qfn- The symmetry-constrained 
variational principle then implies: (i) that the fixed-node 
DMC energy is greater than or equal to the eigenvalue 
of the lowest energy exact eigenstate that transforms ac- 
cording to the same one-dimensional irreducible repre- 
sentation as the trial state; and (ii) that the error in the 
DMC energy is in general second order in the difference 
between the nodal surfaces of the lowest energy exact 
eigenstate that transforms as r r and the trial function. 

If the irreducible representation to which the trial func- 
tion belongs is of dimension greater than one, it is in- 
evitable that some of the symmetry operations from Q 
will change the trial nodal surface. The nodal surfaces 
of the states and ^1 shown in Fig. (|l|) are examples 
of this. The symmetry-constrained variational principle 
need not apply in such cases, because the symmetrized 
state W a used in its derivation need not be zero every- 
where on the trial nodal surface. The delta functions 



produced when the kinetic energy operator is applied to 
the fixed-node pocket ground state 4> a from Eq. (^|) may 
therefore contribute to the expectation value in Eq. JTt|). 

In such weaker version of the symmetry- 

constrained variational principle can be obtained by re- 
analyzing the problem using only the symmetries in the 
subgroup (?fn- The idea is to forget all the symmetry op- 
erations which change the nodes of the trial function, and 
consider only those which leave the trial nodal surface in- 
variant. The symmetry-constrained variational principle 
then applies as long as the symmetries are labeled using 
the irreducible representations of £/fn- 

In summary, a useful DMC variational principle ex- 
ists whenever the trial state transforms according to a 
one-dimensional irreducible representation T r of Q or any 
subgroup of Q . In many DMC simulations, the nodal sur- 
face of the trial state is the same as that of an eigenstate 
of an approximate Hamiltonian such as the local-density- 
functional Hamiltonian. If the state used to define the 
nodal surface is the lowest energy eigenstate with symme- 
try r r of an approximate Hamiltonian with a reasonable 
local potential, the generalized tiling theorem discussed 
in Sec. VI shows that all the nodal pockets are equiva- 
lent by symmetry. The DMC energy is therefore inde- 
pendent of the initial distribution of walkers among the 
nodal pockets. 

The ordinary fixed-node approximation provides a 
good example of these ideas. The many-electron ground 
state is never the overall ground state (which is bosonic) , 
and may be degenerate, in which case we cannot prove 
the existence of a variational principle by analyzing 
the problem using the full symmetry group. We can, 
however, use the permutation group, which is always 
a subgroup of the full symmetry group. The many- 
electron trial function transforms according to the one- 
dimensional antisymmetric irreducible representation of 
this subgroup. The weaker variational principle therefore 
guarantees that the fixed-node DMC energy is greater 
than or equal to the energy of the many-electron ground 
state; and the generalized tiling theorem guarantees that 
the exact many-electron ground state possesses the tiling 
property with respect to permutations. This shows that 
the ground-state versions of the fixed-node variational 
principleEll and tiling theoremH may be regarded as spe- 
cial cases of the more general versions discussed in this 
paper. 

The different members of a set of trial functions form- 
ing a basis for a multi-dimensional irreducible represen- 
tation of Q have different nodal surfaces and need not all 
produce the same fixed-node DMC energy. The strength 
of the weaker variational principle may therefore be opti- 
mized by using specific linear combinations of these basis 
functions. The best linear combinations transform ac- 
cording to one-dimensional irreducible representations of 
subgroups of G, and may be found following the proce- 
dure explained in Sec. M. If this procedure is not carried 
out, trial functions belonging to multi-dimensional irre- 
ducible representations of Q usually have nodal surfaces 



12 



with very little spatial symmetry. In many cases, the only 
symmetry operations that leave the nodal surface invari- 
ant are the elements of the permutation group, and the 
only variational principle that survives is the one relating 
to the many-electron ground state. 
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FIG. 1. The single-electron trial states (a) and (b) 
^2 from Eq. (^). The dashed lines show where the nodal 
surfaces cross the plane of the molecule. The expansion co- 
efficients of $i and ^2 in terms of the tight- binding basis 
functions are also shown. 
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FIG. 2. Angular eigenfunctions Q(8) of the separable 
Hamiltonian of Eq. @ with V(r) = cos(36»)/p 2 . The func- 
tions (a) and (b) are analogous to the tight-binding states 
(a) and (b) shown in Fig. (|l|). The corresponding angular 
eigenvalues Eo are also shown. 
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FIG. 3. Angular eigenfunctions €)(#) of the separable 
Hamiltonian of Eq. @ with V(r) = (cos(30) + 5 cos(66>))/p 2 . 
The functions (a) and (b) are analogous to the tight-binding 
states (a) and (b) shown in Fig. ([l]). The corresponding an- 
gular eigenvalues Eg axe also shown. 



FIG. 5. The total pocket eigenvalues E of th e two nodal 
pockets of the separable example from Sec. 
tonian is that of Eq. (ES) with V(r) 



VIII 



The Hamil- 
cos(3#)/p 2 , and the 
eigenvalues are plotted as functions of the angular half width 
of the nodal pocket centered on 6 = 0. The vertical line shows 
the nodal surface of the exact excited state ^2. 
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FIG. 4. The angular pocket eigenvalues Ee of th e two 
nodal pockets of the separable example from Sec. VIII. The 
Hamiltonian is that of Eq. ( pjjj ) with V(r) = cos(3#)/p , and 
the eigenvalues are plotted as functions of the angular half 
width of the nodal pocket centered on 6 = 0. The vertical 
line shows the nodal surface of the exact excited state ^2- 
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